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ABSTRACT 

When a sufficiently massive satellite (or secondary) black hole is embedded in a gas 
disk around a (primary) supermassive black hole, it can open an empty gap in the disk. 
A gap-opening secondary close to the primary will leave an imprint in the broad com- 
ponent of the Fe Ka emission line, which varies in a unique and predictable manner. If 
the gap persists into the innermost disk, the effect consists of a pair of dips in the broad 
line which ripple blue-ward and red-ward from the line centroid energy respectively, 
as the gap moves closer to the primary. This ripple effect could be unambiguously 
detectable and allow an electromagnetic monitoring of massive black hole mergers as 
they occur. As the mass ratio of the secondary to primary black hole increases to 
q > 0.01, we expect the gap to widen, possibly clearing a central cavity in the inner 
disk, which shows up in the broad Fe Ka line component. If the secondary stalls at 
^ lO^Tg in its in-migration, due to low co-rotating gas mass, a detectable ripple effect 
occurs in the broad line component on the disk viscous timescale as the inner disk 
drains and the outer disk is dammed. If the secondary maintains an accretion disk 
within a central cavity, due to dam bursting or leakage, a periodic 'see-saw' oscillation 
effect is exhibited in the observed line profile. Here we demonstrate the range of ripple 
effect signatures potentially detectable with Astro-H and 1X0/ Athena, and oscillation 
effects potentially detectable with XMM-Newton or LOFT for a wide variety of merger 
and disk conditions, including gap width (or cavity size), disk inclination angle and 
emissivity profile, damming of the accretion flow by the secondary, and a mini-disk 
around the satellite black hole. A systematic study of ripple effects would require a 
telescope effective area substantially larger than that planned for IXO / Athena. Future 
mission planning should take this into account. Observations of the ripple effect and 
periodic oscillations can be used to provide an early warning of gravitational radia- 
tion emission from the AGN. Once gravitational waves consistent with massive black 
hole mergers are detected, an archival search for the FeKa ripple effect or periodic 
oscillations will help in localizing their origin. 

Key words: galaxies: active - (stars:) binaries: close - planets-disc interactions - 
protoplanetary discs - emission: accretion 



1 INTRODUCTION 

Galactic nuclei host supermas sive black holes (> 10® M©) 
iKormendy fc Richstonel Il995l ) and should be the sites of 
mergers of massive black holes. Major and minor galactic 
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mergers in ACDM cosmology should inject supermassive 
and intermediate mass black holes (BH) into the dominant 
galaxy. Following the merger of two BH-harboring galaxies, 
the BHs sink to the bottom of the new galactic potential 
via dynamical friction in approximately a galactic dynam- 
ical timescale (Bcgclman, Blandfo rd fc Recs 198Q). In ad- 
dition to stellar interactions (e.g. jPreto et al.ll201ll ). many 
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studies have shown that gas in the vicinity of the binary 
coul d aid in h ardening the binary down to ^ pc separations 
fe.g.lEscala"et al., 2005i:lDotti et al.ll2007|:lMaver et al.|[2007 : 
Haiman et al.ll2009l: iLodato et al.ll2009l: Icuadra et aliboOE 



Nixon et al.ll201ll : IChapon. Maver fc Tevssieij|201ll '). Black 
holes of intermediate mass can also grow quickly in gas disks 
in active galactic nuclei (AG N) via stellar /stellar remnant 
collis ions and gas accretion (|Levinl l2007l : iMcKernan et al] 
|2012| ) and orbital decay can generate a later merger with 
the central black hole. Massive black hole mergers should 
correspond to some of the strongest sources of gravitational 
waves in the Universe. Since we have yet to detect gravita- 
tional radiation directly, any electromagnetic signature that 
allows us to track massive mergers as they occur will be 
valuable for the study of strong gravity and the growth of 
the largest black holes in the Universe. 

Gravitational torques from satellites in gas disks act 
to repel gas away from the satellite orbit. Sufficiently mas- 
sive satellites c an open empty an nular paps in protoplane- 
tary disks (e.g. iLin & Papaloizou 1986; Pollack et al."l99^; 
IWardlll997l : lArmTtage 2010; Duffell fc Ma cFad ven .2012 . & 
references therein). Analogously, sufficiently massive black 
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tures (jMcKernan et all l2013l ). As the ratio of satellite 
black hole mass to central black hole increases to 10"^ ^ 
q < 1, the gap can widen to form a central cav- 
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therein). A gap-opening black hole in the innermost AGN 
disk, analogous to a gap-opening 'hot Jupiter' in a proto- 
planetary disk, can leave an imprint on a relativistically 
and gravitationally broadened FeKa line profile. The rip- 
ple effects and oscillations that result from a variety of 
black hole merger scenarios: (1) are potentially detectable 
with near-future X-ray detectors such as Astro-H,LOFT and 
IXO/Athena and possibly with present instruments such as 
XMM-Newton , (2) can be used to follow massive merg- 
ers as they happen, (3) provide advance warnin g of p ravi- 
tational radiation from the source l|Kocsis et al. 2008) and 
(4) can test models of extreme gravity, independent of the 
detection of gravitational radiation. However, detecting the 
effects outlined here with present and future proposed mis- 
sions requires serendipity. To systematically survey sites of 
potentially merging massive binaries covering all the effects 
discussed here requires at a minimum, the planned effective 
area of LOFT combined with the planned spectral resolu- 
tion of IXO/Athena. 

In section ij2]we review the broad FeKa line and out- 
line the ripple effect caused by a narrow annular gap around 
the secondary's orbit, deep in the potential well of a super- 
meissive black hole. We also discuss several other different 
plausible configurations for the gas distribution around bi- 
naries. In section §3] we demonstrate the effects of an annu- 
lar gap on the broad FeKa line profile for a range of basic 
observables (angle to observer. X-ray emissivity profile). In 
section §3] we discuss other possible circumbinary gas disk 
configurations for merging BHs including stalled migration, 
damming of the accretion flow, a central cavity and a sec- 



ondary accretion disk around the gap-opening black hole. In 
section iJS]we discuss the gravitational wave signals that may 
correspond to the FeKa broad line profiles. Finally in ^we 
summarize our main results and present our conclusions. 



2 A RIPPLE EFFECT IN THE BROAD FEKa 
LINE 

2.1 The broad FeKa line component 

Broad FeKa lines are observed in several nearby AGN (e.g. 
Nandra et ai] 1 19971 : [Turner et al.ll2002l: iRevnolds fc NowaM 



20031 : 



Braito et al.l l2007l : iMiniutti et all |2007^ . The" nar- 
row core of the FeKa line observed in AGN (e.g. 
lYagoob fc PadmanabhanI bOO'jl') does not correlate signifi 



cantly with outfiowing gas ijMcKernan et al. J2007l ) and so 
presumably originates in fiuorescing cold gas j_far from the 



supermassive black hole l|Shu. Yagoob fc Wand [2OI0I ) . We 
shall not discuss the narrow component further here. How- 
ever, the broad FeKa line component in (AGN) is widely 
believed to originate in fiuorescent Fe deep in the gravita- 
tional potential well of the central supermassive black hole. 
The broad FeKa line profile is complicated by 'horns' due 
to relativistic boosting and a strong red win g due to gravi- 
tational broadening (see e.g. IRevnolds fc Nowak 200 3. for a 
review) . 

There is some recent controversy over the origin 
of the broad component of the line, with some sug- 
gestions that it originates i n complex absorption (e.g. 
I Miller. Turner fc Reeve j |2009| ). However, partial covering 
models do not eliminate the need for a broad FeKa com- 
ponent l|Gallo et al.ll201ll ) and Compton-thick outfiows are 
ruled out in th e innermost regions of the accretion disk 
(|Revnoldsll20l3 ). Furthermore, the ap pearance of the broad 
line during transit events in AGN l|McKernan fc Yagoobl 
119981) is consistent with an orig in in the innermost accretion 
disk (|Weaver fc Yagoobl [l99i ). So it seems we should ex- 
pect a broadened Fe Ka line component originating deep in 
the gravitational potential well around supermassive black 
holes. 

Although the structure of the innermost accretion fiow 
is unknown, it is common to assume a disk structure. 
The innermost parts of a disk mus t lie in the orbital 
plane of a massive merging binary jlvanov et al.l Il999l ). 
In this context, the broad Fe Ka line profile has been 
used t o test models of warped disks jHartnoll fc BlackmanI 
I2OOOI: iFr agile. Miller fc Vandernool |2005|). thick disks 
(IWu fc Wang 2007). accretion rings (ISochora et all [2OI1I ) . 
spiral density waves (jHartnoll fc BlackmanI |2002|) and th e 
presence of a binary with two disks ( Sesana et al.l I2OI2I ) . 
Given the origin of the broad component of the Fe Ka line 
in the innermost regions of the AGN disk, a gap or cavity in 
the inner disk due to the presence of a massive (secondary) 
black hole will have an effect on the broad line profile. We 
discuss this effe ct in detail in t he rem ainder of this paper. 

Following iFabian et al.l (|l989l ). we calculate the 
fiux emitted by a axisymmetric thin disk orbiting a 
Schwarzschild black hole. The ratio of the emitted energies 
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from a point in the disk to the observed energy is given by 
l|Fabian et aLlll989l ) 



(1 + ^) 



-1/2 



1 + 



cos/3 



r(l + tan2f)/M - l]i/2 



(1) 



where (r, (j>) is the location of the emitting point in the disk 
(with r in units of rg = GM/c? , the gravitational radius). In 
eqn. ([T]) M is the black hole mass, P is the angle between the 
disk plane and the plane of the photon trajectory, ^ + 7r/2 
is the angle between the emission of the photon and the 
line connecting the emitting poi nt to the black hol e. The 
quantity cos/3 may be written as (|Fabian et al.lll989h 

cos^sinS 



cos/3 — 



Icos^O + cos^^sin 



11/2 



(2) 



with 6 the angle between the observer and the disk. The 
specific intensity in the frame of the emitting plasma is ap- 
proximated as a delta function /em = e5[E — Eo) where Eo 
is the line rest energy. The structure of both the hot corona 
(source of the X-ray continuum) and the inner disk (source 
of the fluorescing Fe) is unknown, so e the disk emissivity 
is simply assumed to be oc {r/M)~'' . Therefore the specific 
flux measured by the observer from po rtion rdrdij) of the 
disk is then given by (|Fabian et al.lll989l ') 



dFnha — 



do, 



■rdrdcj) 



(3) 



(1-Hz)3^""a(r, cj>) 

where the solid angle dO, = Z?~^bdbd(/>' where D is the dis- 
tance to the disk from the observer, b is the photon impact 
parameter and 

d<i>' cosO 



cos^</> -f sin (Acos^ 



In the weak field limit, the photon geodesies can be ap- 
proximated with straight lines, which means that not all 
relativis tic effects are fully calcu lated at < 20rg in this ap- 
proach (|Beckwith fc Dong|2004l ). By integrating eqn. ^ as 



-^obs — 



dFohs + 



(5) 



where the disk extends from rinnor to router, with a gap exca- 
vated from radii ri to r2, we obtain the observed broad line 
flux from a disk containing a gap. In the limit as ri — >■ r2, 
and if the disk extends to the innermost stable circular or- 
bits (ISCO), rinncr — >■ nsco, wc rccover the observed broad 
line flux from the total disk. In the limit as rinner — >■ '"i we 
recover the observed broad line flux due to the inner disk 
minus a cavity spanning nsco to r2. 

2.2 Gap-opening in disks and merging black holes 

The disk emissivity which produces the broad Fe Ka line 
component depends on both the structure of the X-ray con- 
tinuum emitting corona plus the gas surface density distri- 
bution in the circumbinary disk. The geometry of both the 
corona and inner disk are poorly understood, particularly 
at the late stages of a merger with the central supermassive 
black hole. The main purpose of this paper is to illustrate 
that different possible gas conflgurations produce very dif- 
ferent broad FeKa line proflles. Although our illustrations 
are based on toy models, the main features should be robust, 
and potentially distinguishable in observations with Astro-H 



and IXO/Athena. In this section, we outline possible mod- 
els for the circumbinary gas, including those with a narrow 
annular gap, a fully empty central cavity, or a central cavity 
with possible minidisks inside the cavity. Throughout this 
discussion, we emphasize the very large existing theoretical 
uncertainties. 

BH binaries are commonly believed to be sur- 
rounded by a thin (H/R < 1) circumbinary disk. 
In models where the binary results from a galaxy- 
galaxy merger, the merger is expe cted to deliver the 
nuclear BHs (e .g. Springcl, Di Matte o~fc Hernguistl l2005l : 
[Robertson et a l.' 2006), along with copious amounts of gas 
( Barnes fc Hernq uist.l99Z). to the central regions of the new 
galaxy. The gas cools efficiently, and loses angular momen- 
tum, creating a thin sub-parsec disk in which the pair of BHs 
are e mbedded (e.g. Dotti et a^. 2009; Hopkins & Quataer^ 
I2OIOI : IChapon. Mayer fc Tevssiei]|201ll '). IMBHs with lower 
masses can also grow quickly in situ in the disk around a 
single supermassive black hole, again pro ducing a pair (or 
more) of BHs embedded in a thin disk (|McKernan et all 
2012). 

A BH binary can exchange angular momentum with 
the disk, distorting the disk's density profile, causing the 
secondary BH to migrate inward. In the limit of a very low- 
mass secondary {q = M2/M1 <C 1), the distortions, in the 
form of spiral density waves, remain in the linear regime 
(and the corresponding Typ e I migration is very rapid; e.g. 
iGoldreich fc Tremaind [l980l ) . In this limit, the FeKa line 
would be hardly affected by the presence of the low-mass 
secondary 0. As the mass of the secondary increases, gas is 
gravitationally torqued away from the orbit, lowering the 
co-rotating gas density and a p articularly rapid form of mi- 
g ration (Type HI) may occur (|Masset fc Papaloizoul l200"3l : 
Mc Kernan et al.]|2011bl ). Once the secondary becomes suffi- 
ciently massive it will open an empty annular gap in the disk, 
analogous to gaps known to exist in protoplanet ary disks. 
This happens above the critical mass ratio (e.g. lArmitag^ 
boiof ) of 



\ 8 J 



5/2 



1/2 



(6) 



where H/r is the geometric thickness (aspect ratio) 
of the disk and a is the disk viscosity parameter 
(jShakura fc SunvaevI [l973 ). A black hole with q ^ 10"* 
should be able to open a gap in a fiducial AGN disk. This 
mass range spans massive stars or black holes of > 10^ Mq in 
accretion disks around 10^ Mq supermassive black holes 
l|McKernan et all I2OI2I ). all the way up to supermassive 
black hole binaries that form long after a major galactic 
merger has occurred. The radial extent of the gap would be 
of order a few Hill radii (±i?H). In this case, the secondary 
can be regarded as a particle in the disk, migrating inward 
on the viscous time-scale (so-called Type II migration). 

When the secondary arrives at the radius where the 
local disk mass is too small to absorb the secondary's an- 
gular momentum, its migration must stall (or at least slow 



^ However, if the local gas mass is low enough, or for low enough 
values of H/r, a in the inne rmost disk, even low-ma ss migrators 
expected in the AGN disk (McKcr nan et al.ll20lial . |2012) could 
open a gap or cavity and affect the Fe Ka line as described here. 
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down). An important point is that "standard" thin disks 
have relatively low mass, and this stalling occurs at fairly 
large binary separations. With the exception of very unequal 
masses q < 0.01, the transition takes place well in the outer 
region of the disk, at r >several xlO'^rg. Even for q <^ 1 
binaries, the transition is located at 

retail = 260(Afi/10'MQ)-2/'(g/10-^)'/' (7) 

Schwarzschild radii, which may lie outside the region of in- 
terest for producing broad FcKa lines depending on router 
(the above equation assumes standard disk parameters; 
see eqn. 25 in iHaiman et al.l l2009l for the full expression). 
Once inside this stalling radius, the secondary is no longer 
able to maintain steady Type II migration on the vis- 
cous time-scale, and the evolution of the binary-fdisk sys- 
tem from this point onward becomes very poorly known. 
The inner disk - whose viscous time is short - is be- 
lieved to drain onto the primary, creating a central cav- 
ity. Continued accretion from larger radii causes a pile- 
up of gas just out side the secondary's orbit, analogous to 
a dam in a river (ISver fc Clarke .19951: llvanov et al.l 1 19991: 



Milosavlievic fc PhinnevI l2005l : I Chang et al.ll2010l : iRafikovl 
20121 ). If the dam is 100% efficient in continuing to block 



the gas arriving from large radii, then the binary, at small 
separations, will be surrounded by an empty inner cavity, de- 
void of any gas. Such inner cavities have been seen in many 
numerical simulations (start ing with the seminal work by 
lArtvmowicz &: Lubowl l|l994l )). However, simulations do not 
follow the system on long enough timescales to determine 
whether such a cavity could persist for many viscous times. 
In the face of ongoing accretion, an initial cavity may be 
gradually refilled (if the "dam" is porous), or else the gas 
accumulate d outside the orbit ma y eventually cause the dam 
to "burst" (|Kocsis. Haiman fc Lo cb 2012a,bl). 

An empty cavity would of course strongly modify the 
shape of the FeKa line. However, this picture is over- 
simplified. Many numerical simulations imply that gas 
can enter such a cav ity in narrow coUimatcd streams 
dArtymowicz &: Lubowlil996 : Hav asaki, Mineshige fc Sudou 
20071: iMacFadven fc Milosavlievi3]2008l: ICuadra et al.ll2009l : 



Roedig et al.ll201ll . 120121 : IShi et aLiboTll ). Such streams may 
feed a "minidisk" around the seconda,ry and perhaps also 
the primary BH (|Havasaki et all l2008l : ISesana et 311120121 : 
iFarris et al.ll20l3 ) . A minidisk around the secondary may be 
favored, but the details depend on the angular momentum 
and shocks in the material in the streams. The streamers and 
minidisks could create some additional FeKa emission. For 
exa mple, if both d i sks em it, there would be a double FeKa 
line H esana et all (|2012i )). On the other hand, it is unclear 
whether such streams and minidisks may persist down to 
small bin ary separations. At separations of < 2QQrg<f (see 



eqn. 30 in [Haiman et al.ll2009l for the full expression, includ- 
ing dependencies on disk parameters), the binary BH emits 
GWs efficiently and is driven to merger by this GW emis- 
sion. I t effectively decouples froin the d isk (|Liu. Wu fc Cad 
(|2003t l: iMilosavlievic fc PhinnevI (|2005')) and "runs away". 
The merger timescale is given by (|Pe tcrs 1964) 



TGw ~ 10 yr 



ondary orbital eccentricity. For a secondary on a circular- 
ized orbit with g = 3 x 10"^ located at a = (lO)lOOrg in 
an AGN disk, tgw ~ (30yr)0.3Myr around a 10** M© SMBH 
and TGW ~ (lyr)O.OlMyr around a 10*^ M© SMBH respec- 
tively. Simulations in the relativistic regime found that gas 
streams can follow the binary to small ra dii (~ few Ra\ e.g. 
iNoble et al]|2012l : lFarris et aLll20lj] . l2012l ). but these simula- 
tions start with small binary separations, with gas assumed 
to have followed the BH to these i nitial separations. 

As emphasized recently by iKocsis. Haiman fc Loebl 
(|2012al lbll. even apart from possible streamers and mini- 
disks inside the central cavity, the formation of the cav- 
ity itself, at large radii, is still uncertain. The coupled 
time-dependent processes of the secondary's migration, in- 
ner cavity formation, and damming-up of gas outside the 
secondary's orbit, has not been modeled self-consistently, 
even in one-dimensional calculatio ns. However, using self- 
consiste nt steady-state solutions, iKocsis. Haiman fc Loebl 
(|2012bl ) have argued that in many cases, pile-up can cause 
a dam overflow at large binary separations. In this case, 
when the binary arrives at small radii relevant to the FeKa 
line, the gas configuration would presumably resemble a 
filled disk, with an annular gap (of order the Hill radius). 
Although the above is yet to be demonstrated by a self- 
consistent calculati on, it is at least suggested by recent sim- 
ulations of .Baruteau et al. I (EOI^). This is the only existing 
simulation that looks at small binary separations, in the 
GW-driven regime, and does not manually insert an empty 
cavity in the beginning. They show that the secondary moves 
inward, "ice-skating" across the inner disk (with the inner 
disk material continually crossing the secondary's orbit out- 
ward, on horse-shoe orbits). This simulation is in 2D, and 
assumes a constant kinematic viscosity physics, which may 
not be accurate. Nevertheless, it suggests that if the cavity 
filled in at large radius, it will remain filled until the end, 
except for a narrow annular gap around the secondary. 

We are motivated by the above discussion to consider 
the impact on the Fe Kq line of four different circumbinary 
gas disk geometries, all variations on a standard single-BH 
accretion disk: (i) a narrow annular gap around the sec- 
ondary, (ii) an empty central cavity, (iii) a central cavity 
surrounded by an overdense ring ('dam'), and (iv) a central 
cavity with a small, circumsecondary accretion disk ('mini- 
disk') inside it. We regard (i) as our fiducial case, addressing 
it first in i]2.3l and SjS] We then examine the alternative sce- 
narios (ii)-(iv) in 



2.3 A ripple effect in the broad FeKa line 
component 



Ma 



Ml 



O.OOlpc 



where M2 is the mass of the secondary. Mi is the mass of 
the primary, a is the binary separation and e is the see- 



As a gap-opening black hole migrates inward in an AGN 
disk, the empty gap spans parts of the disk that produce 
ever bluer and redder parts of the broad Fe Ka line. There- 
fore we expect flux 'notches' to be removed from the broad 
line profile. As the annulus moves inward on the disk vis- 
cous timescale, we expect the notches to ripple red-ward 
and b lue-ward respect i vely f rom the line centroid. Fig. [T] 
(from iMcKernan et al.l (|2013h ) demonstrates this ripple ef- 
(l_e2)'^/2(8) feet for a gap-opening black hole located at circularized 
orbit radius R centered in an empty gap of width 2Rh 
where Rh = [q/if^^R. The black hole migrates inward 
with the gap from lOOrg to 6rg. We assumed g = 3 x 10~^ 
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spin (e.R. Dovciak et al, l2004l : [Brenneman fc Revnoldsll2009l : 



Energy 



Figure 1. From [McKernan et alj l|2013h . showing Iron line flux 
in arbitrary units versus Energy in units of keV, binned at ap- 
proximately the energy resolution (~ 7eV) expected for Astro-H. 
The curves demonstrate the progression of the 'ripple effect' in 
the Fe Ko line profile given the in-migration of a gap-opening 
black hole (mass ratio 5 = 3 X 10~^) located at R in the inner 
disk. The empty gap has width 2Rh where Rh = (q/i)^/^R. The 
disk has an assumed X-ray emissivity radial profile of r~^'^ and 
the angle to the observer's sightline is 9 = 60° (where 6 = 0° is 
face-on). The solid black curve indicates the unperturbed Fe K« 
line profile (i.e. no gap). The turquoise solid curve corresponds to 
the same Fe Ka line but with an annulus at 90 it 9rg in the inner 
disk. Dark blue solid curve denotes an annulus at SOitSr^, Green 
solid curve corresponds to an annulus at 20 ± 2rg and the dashed 
red curve corresponds to an annulus at 10 it Ir^. 



for ease of calculation, the disk is oriented at 60° to the 
observer's hne of sight, the X-ray emissivity profile goes 
as r~'^'^ and the broad line flux is assumed to come from 
< lOOrg in the disk (i.e. nnnor = 6ro, router = lOOrp). If we 
change router to larger radii (e.g . iNiedzwiecki fc Zvckl200 
iNiedzwiecki fc Mivakawall201Cl ). the broad hne profiles de- 
picted throughout this paper would appear more centrally 
peaked with 'horns' closer together, resulting in a slight de- 
crease in the prominence of the notches and oscillations 
described here. If router is larger, we become sensitive to 
stalling radii, gaps and cavities further out in the accre- 
tion disk. Exactly how sensitive to detection of the effects 



discussed here depends on the line emissivity (oc r~ ) pro- 
file. If the emissivity profile is steep (k > 3) then most line 
flux originates close to the BH and we have low sensitivity 
to features in the disk much further than ~ lOOrg. If the 
emissivity profile is shallower (k ~ 2), for large router, the 
line 'horns' correspond to material around ~ 100 — ISOr^, 
so gaps at these radii could be detected. In general, our 
simple assumption that the broad line flux comes from 
< 10^ rg), is reasonably consistent with fits to man y observed 
broad line components (jRevnolds fc Nowaklliooi ) . Likewise, 
in the following discussion we do not correctly treat general 
relativistic effects such as light-bending, blurring and BH 



IPauser et al. 2OIOII2OIC ). These effects, if properly treated, 
will tend to smear out sharper notches in the red wing of 
the broad FeKa line and modify line flux originating from 
< 20rg. Even though we ignore the many possible confound- 
ing effects (e.g. continuum absorption, the location and vari- 
ation of the continuum source, broken axi-symmetry), we 
expect that the kinds of spectral features here are broadly 
robust and allow us to take a first step in characterizing dif- 
ferent scenarios of merging massive black holes in the inner 
accretion disk. 

From Fig. [T] the progression of the ripple effect is quite 
clear. As the satellite black hole and gap migrate into the 
inner lOOr^ of the disk, a flux deficit appears (turquoise 
solid curve) compared to the unperturbed profile (black solid 
curve). As migration continues, matching notches are re- 
moved from the broad red wing and the blue wing of the 
profile, with the deepest notches in the blue wing. The twin 
notches move apart, one red-ward, the other blue-ward fur- 
ther from the line centroid (6.40keV) in a distinctive mono- 
tonic progression. As the gap moves inward in the disk from 
SOrg (dark blue curve) to 20rg (green curve) to lOrg (red 
curve), the blue notch moves from 7.1keV to 7.6keV in the 
blue wing, and the red notch moves from 5.6keV to 4.2keV 
in the broad red wing. For a SOOOM© gap-opening IMBH 
around a 10® Af© supermassive black hole, the progression 
from green to red curve in Fig. [1] via GW emission takes 
(eqn. [5]) ~ 16yrs and from red curve to final merger takes 
~ lyr. Because of the shorter timescales, low mass AGN 
exhibiting broad Fe Ka line profiles (e.g. MCG-6-30-15) are 
the prime targets for observations searching for this ripple 
effect, and the massive mergers that cause it. 

AGN that display this ripple effect in the broad Fe Ka 
line will be some of the most luminous sources of gravita- 
tional radiation in the Universe. By searching for the latter 
stages of a progression as in Fig. [l] we can follow the fi- 
nal stages of massive mergers electromagnetically. The final 
stages of a ripple effect in the broad FeKa line will provide 
advance warning of an outburst of gravitational waves from 
this AGN. Conversely, once gravitational waves consistent 
with an IMBH-supermassive black hole merger are actually 
detected, an archival search for the FeKa ripple effect can 
be used to localize the source. 

We expect departures from the predicted progression in 
Fig. \T\ at very small radii either as the disk profile changes 
or as the gravitational radiation luminosity increases dra- 
matically and perturbs the disk (a and H/r may change 
considerably). So the behaviour of the broad component of 
the Fe Ka line may be complicated in the final stages of 
the merger (in the limit of small disk radii). Variations from 
the predicted profile at small radii will include the effects 
of gravitational radiation on the accretion disk. So track- 
ing the ripple effect may allow us to directly test models 
of massive mergers and gravitational radiation, without ac- 
tually detecting gravitational radiation. The Fe Ka band 
can also be spectrally complicated by emission lines. By 
monitoring the monotonic change in energy of the blue and 
red dips, we can break the degeneracy associated with high 
energy Fe emission lines superimpose d on the FeKa com- 
plex (e.g. iMcKernan fc Yaaoobll2004h . Curious 'notch'-like 
features h ave been observed in broad Fe Ka line compo- 
nents (e.g. lYagoob fc Serlemitso"3l2005h . but the collecting 
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area and the energy resolution of the high energy trans- 
mission gratings on board Chandra are insufficient to re- 
solve these features and follow them over time. The broad 
component of t he Fe Kg line can vary quite significantly 
over time (e.g . Iwasawa et al.l Il996l : IdeMarco et al.1 12OO9I : 
ISvoboda et alll2012l ). but it has been difficult to determine 
the nature of the variability given detector limitations. 

The energy resolution of the micro calorimeter planned 
for the Astro-H missio n will be ~ 7eV in the FeKa band 
l|Takahashi et al.l I2OI0I '). This exceptional resolution com- 
bined with moderate collecting area should be sufficient to 
distinguish between absorption features and narrow, mono- 
tonically migrating notches of the kind depicted in Fig. [1] 
in lo ng, repeated observa tions of nearby. X-ray bright AGN 
(e.g. iNandra et a"i]|l997l ). However, it is important to note 
that Fig. [1] displays only the predicted broad line fiux. To 
systematically obtain broad line spectra that can distinguish 
between the different cases in Fig. [T] requires a minimum 
telescope effective area of that planned for LOFT (~ lOm^ 
at 6.4keV), together with the planned spectral resolution 
of (~ 5 - 7eV at 6 4key ) for Astro-H and 1X0/ Athena. 
IVaughan fc FabianI (|2004l ') obtain ~ 18,000 counts in the 
broad red wing of the Fe Ka line in a ~ 320ks observation 
of MCG-6-30-15 with the EPIC detector on board XMM- 
Newton , yielding ~ 600 average counts at the limiting en- 
ergy resolution of ~ O.lkeV per bin and very high signal to 
noise. An observation with a future Large Area High Reso- 
lution (L.A.H.R.) telescope with the effective area of LOFT 
would get ^ 6,000 average counts per ~ O.lkeV bin at very 
high signal to noise. In this case, observations with LAHR 
would have the signal-to-noise of XMM-Newton EPIC at 
binnings of ~ lOeV, which would be sufficient to at least 
systematically track the monotonic changes in the energy of 
the notches in Fig. [1] In the remainder of this paper we plot 
all figures at the ~ 7eV energy resolution expected for Astro- 
H (and 1X0/ Athena), to illustrate the details that might be 
observed in the Fe K broad line complex at sufficiently high 
signal-to-noise with future detectors. As we shall show be- 
low, the EPIC-pn detector on board XMM-Newton is most 
useful in testing models of oscillations in the broad Fe Ka 
line due to secondary mini-disks, until the advent of LOFT. 
In the case of oscillations, high energy resolution is less im- 
portant than binned-up count-rates and the superior effec- 
tive area of EPIC makes this the present detector of choice 
and we recommend that observers search for serendipitous 
oscillations in the Fe K band when observing with this in- 
strument. The planned 1X0/ Athena mission will have both 
high energy resolution and large effective area. However, sys- 
tematic testing of models of massive black hole mergers via 
the broad Fe Ka line will require Large Area High Resolu- 
tion {LAHR) telescopes with minimum effective area that 
proposed for LOFT, together with the spectral resolution 
planned for 1X0/ Athena. We recommend that future mis- 
sion planning take this into consideration. 

A related issue to detectability is the expected rate of 
occurence of binaries (of mass ratio q) in AGN. There are at 
present a small number of bright, nearby AGN that exhibit 
broad Fe Ka lines. Even with an order of magnitude increase 
in telescope effective area, our sample size would be ~ 10^ 
objects. The probability of seeing a close binary per AGN 
is N2ticsNohs /tdisk where N2 is the number of secondaries 
(for a given range of mass ratios) in the AGN disk over its 



lifetime (tdisk), ires is the residence time of the secondary 
at the range of disk radii where we can observe the effects 
discussed here and A'obs is the number of minimal S/N ex- 
posures of the AGN of duration (tmin) in which the effect 
could be detected. For the lowest mass (Type I) migrators 
we might expect N2 ~ lO'' in an AGN disk that lasts for 
tdisk ~10Myrs around a 10* Mq SMBH, where tros ~ lyr 
(|McKernan et aLll2012f ). For a fiducial LAHR telescope, we 
might expect tmin ~ lOks so in a total exposure over a mis- 
sion lifetime of IMs, we would have a 10% chance of de- 
tecting a low-mass migrator in any one AGN. For SMBH 
binaries, t^cs ~ lO^'yr or more, so as long as A'2 for massive 
SMBH is not vanishingly small, we will also have a reason- 
able likelihood of detecting SMBH binaries with future mis- 
sions. F or SMBH-IMBH binari es, N2 is likely to be of order 
several (iMcKernan et al.1 1201^ 1. tres is likely to be ~ 10'*yr 
in the example above, so in a large sample of X-ray bright 
AGN, we may detect a small number of these binaries via 
effects in the broad Fe Ka line. 



3 BROAD FEKa LINE PROFILES: BASIC 
OBSERVABLES 

In this section we discuss the effects of a few basic observ- 
ables on the expected broad Fe Ka line component profile, 
including angle to the observers' sightline(&) , assumed X- 
ray emissivity profile (oc r~^) and the width (Ar = r2 — ri) 
of the annulus gap. In practice, the broad profiles discussed 
here will be those residuals in the Fe K band after the sub- 
traction of a narrow (very slowly varying) Fe Ka compo- 
nent due to cold fluor escing material far from the supe rmas- 
sive black hole (e.g. lYagoob fc PadmanabhanI l2004f ) . The 
Fe K line complex actually includes two components (Kai 
and Ka2 separated by 13eV), centered at 6.40keV and an 
additional Fe K/3 component at 7.06keV. The Fe K/3 line 
should have ^ 13% of the normalized flux of the FeKa line 
(|Palmeri et al.l |2003|) and so adds a small percentage change 
in the proflle of the overall broad Fe K line, particularly on 
the blue-wing. However, the line profiles that we generate 
here, distinguishing between various merger situations, are 
not materially effected by the addition of a broad Fe K/3 
component, so we shall not discuss that component further 
here. Furthermore, although individual highly ionized line 
components can appear in Fe K band spectra, these are suf- 
ficiently narrow that we can hope to disentangle their signa- 
ture from broad components with Astro-H. For the purposes 
of the present discussion, in order to understand the sense in 
which different binary disk geometries affect the broad Fe K 
line component, we shall ignore the contribution of highly 
ionized Fe to the broad Fe Ka line complex. 

We calculate the broad FeKa line profile around a 
Schwarzschild black hole using a modified version of the 
algorithm that ge nerates the XSPEC diskline model 
(jPabian et al.l [l989l ), described in i)2.1l above. We modify 
the DISKLINE algorithm by subtracting annuli of particular 
widths {At = r2 — ri) or by introducing additional emis- 
sion due to damming of the accretion flow or accretion disks 
around the secondary black hole (see discussion below). For 
ease of comparison of basic observables, in this section we 
discuss an empty gap opened at 20 ± 2rg in an accretion disk 
around a supermassive black hole (the green solid curve in 
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Figure 2. Taking tlie 'ripple' due to an annulus at 20±2rg (the 
green solid curve in Fig. [TJ and clianging the viewing angle (6) 
to the observer. The black and green curves correspond to the 
black and green curves in Fig. [T] respectively (6 = 60° , no gap &; 
gap). The dark blue, red and light blue curves correspond to the 
green curve observed a.t 6 = 20, 40, 80° respectively. All curves are 
binned at approximately the energy resolution (~ 7cV) expected 
for Astro-H. 
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Figure 3. Taking the 'ripple' due to an annulus at 20±2rg (the 
green solid curve in Fig. [TJ and changing the X-ray emissivity 
profile (r~*). The black and green curves correspond to the black 
and green curves in Fig. [T] respectively (with k = 2.5). The dark 
blue, red and light blue curves correspond to A; = f .8, 2.0, 2.8 
respectively. All curves are binned at approximately the energy 
resolution 7eV) expected for Astro-H. 



Fig- HI: where all the broad line emission is assumed to come 
from disk radii spanning [rinncr ~ 6rg, router ~ lOOrg]. We 
then change basic parameters in order to understand the 
sense in which the ripple effect changes. In section i|3.1l we 
consider the impact of disk inclination angle (9) on the ob- 
served ripple effect. In section i]3.2l we discuss how changing 



X-ray emissivity (oc r~ ) as a function of radius as well as 
empty gaps of different widths (Ar) alters the ripple effect. 
In section i]3.3l we briefly discuss how the spin of the central 
black hole might change observations of the ripple effect. 



3.1 Changing inclination angle 

There is a relatively wide range of viewing angles (6) possible 
for observations of Scyfcrt 1 AGN nuclei that display broad 
FeKa lines (|Nandra et al.lll997l ). Indeed if disks are warped 
close to a spinning black hole, even f ace-on AGN may con- 
tain highly inclined inner disks (e.g. Bardeen fc PettersonI 
ll975l : lKing fc Pringle|[200^ : iNixon &: KinEll20ld ). The most 
famous example of a broad Fe Ka line in an AGN (MCG- 
6-30-15) is best-fit with d isk-inclination angle of 9 — 27° 
l|Revnolds fc NowakI 1200^ ). where 6 — 0° corresponds to 
face-on. 

Fig. [2] shows the effect of varying inclination angle (9) 
on the ripple due to the annulus at 20 ± 2rg (green curve) in 
Fig. [1] From Fig. [2] as the angle to the observers' sightline 
is increased to ~ 80° (hght blue curve), the FeKa line width 
increases and the 'notches' become more prominent. The 
'blue' notch shifts away from the H-like Fe xxvi Lya line 
at 6.97keV which prevents possible line confusion. However, 



at this inclination angle, we should expect that intervening 
material (the dusty torus, an outflowing wind or obscuring 
clouds) should absorb softer X-rays, influencing the shape 
of the continuum and therefore the broad component of the 
Fe Ka line. Thus, unless the AGN is 'naked' (i.e. without ob- 
scuring torus or clouds), or the inner disk is highly warped, 
a broad FeKa line at this inclination should be difficult to 
disentangle from absorption effects. 

From Fig. [5] as 9 decreases to nearer face-on (dark blue 
and red curves), the likelihood of obscuring structures de- 
creases, but the blue wing of the Fe Ka shrinks back towards 
the rest-frame line energy (6.40keV). As a result the 'blue' 
notch becomes much narrower (and harder to detect with 
Astro-H). The 'red' notch remains prominent and should 
stand out clearly in observations with Astro-H, until obscur- 
ing material starts to absorb the continuum significantly at 
this energy. 

3.2 Changing X-ray emissivity & gap width 

The geometric structure of the corona (source of the X-ray 
continuum) and that of the inner disk (source of the fiuores- 
cent Fe) are unknown. Therefore the X-ray emission from 
the inner disk is simply parameterized as a radial power- 
law form (?■"''), w here likely ranges for the emissivity span 
fc ~ 1.5 - 3 (e.g. [Reynolds fc Nowak 2003j . The most fa- 
mous example of a broad Fe Ka line is best-fit with X-ray 
emissivity profile of fc ~ 3 l|Revnolds fc NowakI |2003| ). 

Fig. O shows the effect of varying the index (k) of the 
X-ray emissivity powerlaw ((r/M)"*"), on the ripple due 
to the annulus at 20 ± 2rg (green curve) in Fig. [T] From 
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Figure 4. Taking tlie 'ripple' due to an annulus at 20 ± 2rg (the 
green solid curve in Fig. [TJ and changing the annulus width. The 
black and green curves correspond to the black and green curves 
in Fig.[T]respectively (with annulus located at 20±2rg). The red, 
dark blue and light blue curves correspond to annuli located at 
20 ± Irg , it4rg , it6rg respectively. All curves are binned at ap- 
proximately the energy resolution (~ 7eV) expected for Astro-H. 



around a maximally spinning black hole, the Fe Kg line will 
extend rcdward to very low energies (^ 3keV|l wasawa et al.l 
(Ji}96)). In the final stages of merger, a gap-opening black 
hole could leave an imprint (at < 6rg) on a very broad Fe Ka 
line. Around a low mass, rapi dly spinning black hole s uch as 
in MCG-6-30-15 (~ IO^Mq, iMcKernan et al.1 (|2010D '). the 
very late stages of a ripple effect could be observed, but 
would last only a few weeks. A gravitational wave trigger 
may therefore be useful in finding these late stage events. 
The possibility is intriguing, since tracking a massive merger 
electromagnetically to the very final stages would be a pro- 
foundly important advance in this field. 



4 TESTING MERGER MODELS WITH THE 
BROAD FEKq LINE PROFILE 

In this section we examine alternative models for the ge- 
ometry of the innermost AGN disk, as discussed in § 12.21 
above, and ask how we might distinguish between these dif- 
ferent models using the broad Fe Kq line component. In 
>i4.1l we discuss cavity formation and how we might test the 
formation and presence of cavities using the broad FeKa 
line profile. In ^4.2l we discuss observational consequences of 
'damming' the disk due to pile-up of infiowing gas behind a 
stalled, migrating secondary black hole. In i]4.3l we discuss 
the possibility of a disk around the secondary black hole and 
the efi'ect this has on the overall broad Fe Kq line profile. 
In reality, a combination of all three effects is likely to be 
present simultaneously. 



Fig. |3l as fc gets larger, more and more line flux comes from 
the innermost regions of the disk, the larger the notches in 
the line become, and the more pronounced the ripple ef- 
fect. If fc > 3 due to the special case of magnetic torquing 
l|Revnolds fc Nowakll20o3 ). the late-stage ripple efi'ect will 
become even more apparent in observations. 

Fig. |4] shows the effect of varying gap width on the rip- 
ple due to the annulus at 20 ± 2rg (green curve) in Fig. [T] 
From Fig. 31 as we should expect, the notches in the Fe Ka 
line profile become more pronounced (detectable) as the gap 
width increases. In ^ below we discuss in more detail the 
effect on the broad Fe a line profile of increasing the gap 
width via disk drainage to form a cavity, damming the in- 
flowing gas, and the possibility of an accretion disk around 
the secondary black hole. 

3.3 Changing black hole spin 

In the discussion above, we have discussed broad Fe Ka line 
profiles that originate in fluorescent Iron deep in the poten- 
tial well of a non-spinning (a = 0, Schwarzschild) black hole. 
In this case, the accretion disk extends inward as far as the 
innermost stable circular orbi ts (ISCO), which for a = 
corre sponds to nsco ~ 6rg l|Bardeen. Press Teukolskvl 
I1972D . Maximally prograde spinning black holes (a ~ M) al- 
low the inner edge of the accretion disk to extend practically 
to the event horizon risco ~ rg. By contrast, black holes 
with retrograde spin compared to the disk (a ~ —M) have 
risco ~ Qrg. For all except the final stages of merger, the 
spin will make little difference to the ripple effect. However, 



4.1 Detecting cavities in the inner AGN disk 

Fig. [5] shows the effect of a cavity in the innermost disk 
on the ripple due to the annulus at 20 ± 2rg (green curve) 
in Fig. [T] In Fig. [S] we demonstrate the effect of removing 
(keeping) all the fluorescing Fe between 6 — 18rg in the red 
curve (green curve). Note that while this situation corre- 
sponds to a stalled migrator at 20r-g, GW emission is likely 
to rapidly re-start the stalled binary. We choose this value 
of the separation to simply demonstrate differences in line 
profiles. The effects are apparent on the broad Fe Ka line 
profile, where the red and blue wings in the green curve, out- 
side the notches (at ~ 5.0, 7.4keV respectively), vanish from 
the red curve. One problem that becomes apparent from 
Fig. [5] is that it will be more difficult to detect the effect 
on the Fe Ka line of black holes that create inner cavities 
than gap-opening black holes. This is not only because the 
obvious features (the notches) are gone, but also because 
spectral model fltting involves independent flux normaliza- 
tions for both the continuum and the Fe Ka line. Thus, an 
observation of a line profile best fit by the red curve in Fig.[5l 
could be produced by a satellite black hole in a cavity, but 
could also be produced by a complete disk observed at a 
different angle {6) or with a shallower emissivity profile (e.g. 
compare with the dark-blue curve in Fig. [3] corresponding 
to r~^'*). In this case, model degeneracy can be broken by 
tracking the width of the line profile over time. If a line pro- 
file similar to the red curve in Fig. [5] was observed in an 
AGN, a cavity should decrease in radius and therefore the 
fluorescent line width should increase monotonically on the 
disk viscoust imescale. Measurements from other broad lines 



© 2013 RAS, MNRAS OOO.nifTSi 



Ripple effect in broad FeKa lines 9 



CD 




3 4 5 6 7 B 



Energy (keV) 

Figure 5. Comparing the 'ripple' profile due to an annulus at 
20 ± 2rg (the green solid curve in Fig. [l} with the line profile 
due to a disk minus an empty cavity at r < 22rg. The black and 
green solid curves correspond to the black and green solid curves 
in Fig. [T] respectively. The red curve corresponds to the FeKa 
profile from a disk with the same emissivity (r"*) and inclination 
angle {6), but with an empty cavity at radii < 22rg. All curves are 
binned at approximately the energy resolution (~ 7eV) expected 
for Astro-H. 

in the same AGN may also constrain the angle (6) to the 
observer of the innermost disk and rule out e.g. near face-on 
disks. 

In Fig. [6l we follow the draining of the inner disk due to 
a black hole migrator that has stalled at 50rg. The satellite 
black hole has opened a gap in the disk at 50 ± 5rg We 
assume that the outer gap edge remains fixed at 55rg and 
that the inner disk drains on the viscous disk timescale. In 
the time it takes the inner disk to drain, we assume that 
no gas enters the gap and for now we neglect the effect of 
'pile-up' of gas mass at the outer edg e of the gap (5 5rg). 
The viscous disk timescale is given by (|Armitagdl2O10l ) 



where w is the Keplerian angular frequency. For a puffed- 
up inner disk with H/r ~ 0.1 and a ~ 0.01, the inner disk 
will drain in ~ lO^ra where rd is the dynamical timescale 
at the gap inner edge (in this case 45rg). Note that the 
viscous timescale is shortest at small radii, so we expect 
the inner disk to drain 'inside out', i.e. nnncr — > ri over 
time. Therefore we expect the surface density to first drop 
dramatically at small radii as material gets accreted. The 
surface density of gas in annuli at large radii in the inner 
disk will drop as it is 'smeared out' across interior radii as it 
accretes, but we shall ignore the mean emission from interior 
to Tinner and simply approximate the inside out drainage as 
a migration outward of rinncr to ri . 

Fig. [6] shows the effect as the inner disk drains in this 
simple 'inside out' manner. The black curve shows the in- 



ner disk as for the dark blue curve in Fig. [T] i.e. anrmlus 
spanning 45-55rg, inner disk radius at Gvg. The progression 
in Fig. [S] shows the effect of allowing the inner disk edge 
(rinncr) to migrate outward on the disk viscous timescale. For 
a stalled q = 3x 10~^ secondary black hole at 50rg around a 
W^Mq (10® Mq) supermassive black hole, the inner disk will 
drain on the approximate timescale ~ 0.5(50)yrs, far shorter 
than the timescale (row) for merger by gravitational wave 
emission. The progression from black curve to red curve 
(nnncr ~ lOrg), to grccn curve (nnncr = 20rg) occurs very 
quickly. The disk drainage progression in Fig.|5]spends most 
of its time between the dark blue and light blue curves. Ob- 
servationally, the key characteristics to indicate this 'inside 
out' disk drainage are that the positions of the notches in 
the broad line (due to the gap) remain fixed, but the red and 
blue wings are suppressed monotonically as the line inten- 
sity diminishes. In this case, multiple repeated observations 
using the large effective area of XMM-Newton EPIC could 
be used to collect flux from the red and blue wings of the 
broad line. Even with the limited EPIC energy resolution 
(^ O.lkeV), the equivalent width of the line wings (as well 
as the notch positions) can be measured down to fractions of 
individual exposures (~ lOks intervals). However, we expect 
the monotonic decrease in the equivalent width of the wings 
to be most signiflcant over longer intervals (~ yr) while the 
positions of the notches remains fixed. It is possible that 
disk drainage could occur on timescales faster than fiducial 
estimates of Tq,. For example, the sharp gradients at the 
cavity edge can cause viscous diffusion to be significantly 
faster than t he local viscous timescale for uniform (non- 
piled -up) gas (iMilosavlievic fc Phinnevl[2005l : lHaiman" et al.l 
I2OO9I ). Furthermore, the effective viscosity a has been found 
to increase locally near the cavity edge by more tha n an or- 
der of magnitude (|Shi et al.ll201li : lNoble et al.ll2012l '): a sim- 
ilar increase i s seen in the inner regions of single-BH disks 
(|Penna et all 120121 '). So, even for secondaries around very 
massive primaries, changes in the Fe Kq line due to disk 
drainage could be detectable over a period of several years. 
The reverse proces s of gap-filling (Tanaka & Menou 201^; 
iTanaka et allbOlOl ). starting after decoupling from the disk 
and lasting beyond the MBH merger, can produce a time- 
reversed sequence lasting--^ lOyrs for a IO^Mq SMBH. 

If instead, we assume the drainage is 'outside in', the 
disk depletes first at larger radii and then at smaller radii in 
to risco, to form the cavity, the change in the line profile 
will be different. Outside-in drainage is unlikely, since the 
viscous timescale is shortest for small radii, although a{r) 
and H/r in eqn.|^may change significantly with decreasing 
radius. Fig. [7] shows the effect of outside-in drainage, where 
ri decreases to rinncr. For most of this time, the line profile 
will be somewhere between the black and green curves. The 
final stages of disk drainage will happen very fast and the 
rate of change in the broad FeKa line profile will be large 
(going from red curve to purple curve) . As the inner gap edge 
decreases to 40rg (dark blue curve) , 30rg (green curve) , 20rg 
(red curve), lOrg (light blue curve) and r-g (i.e. empty cavity; 
purple curve). 

If the gap-opening black hole stalls further out, say at ^ 
lOOrg or greater as might be expected (|Haiman et al.ll2009l ). 
the inner disk will drain on the viscous timescale of the inner 
edge of the gap at ^ lOOrg, while the outer disk is dammed. 
In this case, for a g = 3 x 10~^ gap-opening black hole 
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Figure 6. Comparing the 'ripple' profile due to an annulus at 
50 it 5rg (the dark blue solid curve in Fig.[T]l with the line profile 
due to a steadily draining inner disk which drains 'inside out'. 
The black solid curve corresponds to the dark blue solid curve in 
Fig. [T] To mimic an 'inside out' draining inner disk, we held both 
the inner and outer gap edges (at 45, 55rg respectively) fixed, 
while allowing the inner disk edge (rinnc^ = 6rg for the black 
curve) to migrate outward (on the disk viscous timescale) to lOr'g 
(red) , 20rg (green) , SOrg (dark blue) , 40rg (light blue) . For a fidu- 
cial disk with a = 0.01, (H/r) = 0.1 around a 10*^(10*^/0) su- 
permassive black hole, we expect the progression of profiles from 
solid black curve to light blue curve depicted here to last approx- 
imately ~ 0.5(50) years. All curves are binned at approximately 
the energy resolution (~ 7eV) expected for Astro-H. 
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Figure 7. Comparing the 'ripple' profile due to an annulus at 
50 ± 5rg (the dark blue solid curve in Fig. [ij with the line pro- 
file due to a steadily draining inner disk. The black solid curve 
corresponds to the dark blue solid curve in Fig. [T] To replicate 
an 'outside-in' draining inner disk, we held the outer gap edge 
(at 55rg) fixed, while allowing the inner gap edge (ri) to migrate 
inward (on the disk viscous timescale) to 40rg (dark blue), 30rg 
(green), 20r'g (red), lOrg (light blue) and Org (purple). This form 
of drainage is unlikely as the viscous disk timescale is shortest at 
small radii, although a and H/r may change significantly as radii 
decrease. For a fiducial disk with a = 0.01, (H/r) = 0.1 around a 
10''(10*Mq) supermassive black hole, we expect the progression 
of profiles from solid black curve to purple curve depicted here to 
last approximately ~ 0.5(50) years. All curves are binned at ap- 
proximately the energy resolution (~ 7eV) expected for Astro-H. 



stalled such that ri = lOOrg in a a = 0.01, H/r = O.l gas 
disk around a 10® Mq(1O*M0) supermassive black hole, the 
inner disk will drain in ~ 1.5(150)yrs ^ tgw- If the disk 
drains 'inside out' (rinnor ri over time), as we expect since 
the viscous timescale is shortest at small radii, the broad 
line component over that time will follow the progression 
in Fig. |8] assuming that disk gas at ^ lOOrg is held up by 
the 'dam wall'. In our example of a q = 3 x 10~^ secondary 
around a 10®(10**)A'/q primary, the line profile will change 
from the black to green curve over 0.5(50)years and from 
green curve to light-blue curve in l(100)years. If the inner 
disk is draining 'inside out' in this way, the key observational 
characteristic to search for is the increasing suppression of 
the red and blue wings of the line over time. This model can 
be tested using archival and ongoing observations with the 
large effective area of XMM-Newton EPIC. In this case, the 
binned-up wings of the broad component will decrease in 
a consistent monotonic pattern over time (quickly at first, 
then more slowly). If instead, the a and H/r profiles of the 
disk mean that inside-out drainage does not occur, and an 
'outside-in' disk drainage somehow manages to happen, the 
red and blue 'horns' of the line profile will ripple outward 
from the line centroid energy as the overall broad component 
intensity diminishes. In this case, the low energy resolution 
of XMM-Newton can capture the migration of the horns 



red-ward and blue-ward, but only over long timescales. The 
higher energy resolution of Astro-H and IXO / Athena could 
resolve horn migration over much shorter timescales. 

Note that if the X-ray continuum originates in a hot 
corona concentrated mostly above the innermost disk then 
during and after cavity formation, we should expect some 
fraction of the overall X-ray continuum to decline. If the 
inner disk drains 'inside out' to form a cavity, we should 
expect a rapid drop of the X-ray continuum together with 
a strong decrease in broad line intensity (most pronounced 
in the wings). Observing such an effect would allow us not 
only to test models of massive mergers, but also allow us to 
constrain the fraction of the X-ray continuum that originates 
in the central ~ lOOrg . If the disk drains 'outside in' (much 
less likely), we expect the drop in the continuum to match 
the rippling apart of the horns of the Fe Ka line. 

4.2 Detecting pile-up in AGN disks 

In the discussion of disk drainage above, we ignored the 
effect of the continuing infiow of gas from the outer 
disk. As the migrating secondary stalls and the disk in- 
terior to it drains on the viscous timescale, mass is still 
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Figure 8. The change in the broad FcKa line profile due to the 
draining of the entire inner disk due to a stalled gap-opening 
satellite black hole, whore ri = lOOrg . To replicate a draining in- 
ner disk, we assumed the broad FeKa emission originates entirely 
from within lOOrg, with 9 = 60° and r~^'^ X-ray emissivity pro- 
file as assumed above. We do not see 'notches' as in Fig. |6] since 
the gap does not imprint itself on the inner disk < lOOrg. The 
solid black curve corresponds to the broad FeKa profile from the 
inner disk extending from lOOrg inward to 6rg. Since we expect 
the disk viscous timescale to increase with radius, we assume the 
disk drains in an 'inside out' manner, so ringer — > ti s.nd we 
ignore emission from gas 'smeared out' within Tinner- As the in- 
ner edge of the disk increases radially to 20rg (red curve), 40rg 
(green curve), 60rg (dark blue curve), SOvg (light blue curve), 
the broad component of the line decreases in magnitude and the 
blue and red wings are increasingly suppressed. For a fiducial 
disk with a = O.Ol,{H/r) = 0.1 around a 10^(10* Mq) super- 
massive black hole, we expect the progression of profiles from 
solid black curve to light blue curve depicted here to last approx- 
imately ~ 1.5(150)yrs. Two-thirds of this time is taken up with 
the progression from green curve to light-blue curve. All curves 
are binned at approximately the energy resolution (~ 7eV) ex- 
pected for Astro- H. 

flowing inward in the disk. Gas should build up at the 
outer gap or cavi ty edge on the visco us timescal e (tq) at 
that disk radius (ISver'fc Clarkd 1 19951 : iRranov et aP 1 19991 : 
iKocsis. Haiman &: Loebll2012al ). We can think of the outer 
gap radius or cavity edge as a dam holding back the inflow 
of gas. Here we discuss how damming inflowing gas at the 
outer wall of the gap or cavity has an effect on the observed 
broad Fe Ka line component. 

We created a naive toy model to represent the damming 
of inflowing gas at the outer gap wall and we considered 
two simple cases. First, as the inner disk drains onto the 
primary, creating a cavity in the inner disk (as in Figs. [7] 
and [8] in ^4.1|) . gas piles up at the outer gap wall. Sec- 
ond, gas piles up at the edge of a pre-existing cavity. In 
both cases our toy model assumes a simple, uniform en- 
hancement of disk emission in a annulus at the disk trun- 
cation radius. Pile-up behind the dam is unlikely to be uni- 
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Figure 9. The change in the broad FeKa line profile due to 
draining of the entire inner disk ('outside-in') due to a stalled 
migrating satellite black hole at 50rg , while gas piles up at the gap 
outer edge. The solid black curve corresponds to the broad FeKa 
profile due to the inner disk with a cavity spanning 45 — 55rg 
(same as black curve in Fig. O. The red curve corresponds to 
ri = 30rg , r2 = 55rg with an enhancement of emission by X 2 
in the annulus spanning 55 — 60rg . The green curve corresponds 
to ri = 10rg,r2 = 55rg with an enhancement of emission by 
X 4 in the annulus spanning 55 — 60rg . All curves are binned at 
approximately the energy resolution (~ 7eV) expected for Astro- 
H. 

form in reality, particularly for large values of q. A dense 
tidal bulge tracking the secondary can form on the dam 
wall (Shi ct al." 2011; Roodis ot al. 201l|; iNoble et al.ll2012l : 
iD'Oraz io, Hai man fc MacFadyen 2012|), which will orbit at 
the orbital time for the cavity wall. This can lead to an os- 
cillation in the effect discussed below but we shall return to 
this in future work. In order to move the stalled secondary, 
we should expect a pile-up comparable to the mass of the 
secondary. We translated this into a simple density enhance- 
ment over a standard thin disk of a factor ~ 2 — 5 distributed 
uniformly within ^ 10% of the cavity edge. This is roughly 
cons istent with profiles of the pile-up seen in simu lations 
(e.g. iMacFadven fc Milosavlievi3l2008l : ICuadra et aLll2009l ) , 
although over time we might exp ect a S(r) ex build 
up to large radii behind the dam (|Kocsis. Haiman fc Loebl 
l2012al ). We shall investigate the effects of different types of 
'damming' in future work. Here we are not attempting to 
reproduce details of the damming; we neglect disk heating, 
pressure and expansion. Rather, we simply want to under- 
stand the sense in which the broad Fe Ka line component 
changes when we enhance the Fe Ka emission at the disk 
edge in a naive manner. 

In the first case, where the inner disk at < ri is drain- 
ing onto the primary as in Fig. [T] we assume that there is 
an enhancement of the emission from gas in an annulus im- 
mediately outside r2. Fig. [9] shows the inner disk draining 
'outside in', but now we uniformly enhance Fe Ka emis- 
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Figure 10. The change in the broad FeKo hne profile due an 
inner disk truncated by a cavity at ri = r2 = 55rg, while gas 
piles up at the gap outer edge. The solid black curve corresponds 
to the broad FeKo profile due to the inner disk truncated by a 
cavity at < 55rg (same as purple curve in Fig. . The red curve 
corresponds to an enhancement of emission by X 2 in the annulus 
spanning 55 — 60rg . The green curve corresponds to enhancement 
of emission by X 5 in the annulus spanning 55— GOrg . All curves are 
binned at approximately the energy resolution (~ 7eV) expected 
for Astro-H. 

sion in an annulus spanning 55 — 60rg immediately outside 
r2 to model pile-up. As the inner gap edge (ri) decreases 
to 30rg (red curve), we add a uniform (arbitrary) x2 en- 
hancement of the Fe Ka emission from disk radii spanning 
55 — 60rg. As the inner gap edge decreases to ri = lOr-g 
(green curve), the Fe Ka emission from disk radii spanning 
55 — GOvg is enhanced by a uniform (arbitrary) x4 factor. 
By comparing Fig. |9] with Fig. [71 we can see that substan- 
tial damming of the outer disk has a noticeable effect on 
the broad component of the Fe Ka line. In particular we 
see that the 'horns' of the broad component each become 
double-pronged as the 'pile-up' beyond increases. Obser- 
vationally, we should expect the double-pronged horns to 
persist until the dam bursts. If the dam burst is sudden and 
catastrophic, we expect the extra horns in the line profile 
will ripple outward to redder and bluer energies respectively 
as the 'overdense' material collapses inward. If instead, there 
is a slow leak from the dam, the flux in the extra horns in the 
line profile will decrease and become dispersed to energies 
redward and blueward of the horns. 

In the second case, we assume that a cavity has been 
excavated, truncating the disk at ri — r2 = 55rg. We also 
assume the outer disk is damming inflowing gas in an annu- 
lus spanning 55 — 60r-g . Fig. [TO] shows the eflect over time 
of uniform pile-up at the edge of the cavity (r^). The black 
curve shows the broad line from the region outside the cavity 
(> 55rg), with no damming of the accretion flow (equivalent 
to the purple curve in Fig. [7]). The red curve shows an en- 
hancement of x2 over 'normal' emission in an annulus span- 



ning 55 — 60rg. The green curve shows an enhancement of 
X 5 over 'normal' emission in the annulus at 55 — 60r-g . From 
Fig. IIUI the key observational characteristic of a 'damming' 
of the accretion flow is a double-horned broad Fe Ka line 
component (as in Fig. [5)). This is because the broadest com- 
ponent in the observable Fe K line complex is emitting sub- 
stantially more than the surrounding disk. In this case, the 
blue- and red-shifted horns of the line from the innermost 
annulus stand out relative to the blue- and red-shifted horns 
of the line from annuli further out in the disk. At some point, 
when the pressure is large enough, the dam will overflow. If 
the overflow is continuous, we expect a dynamic equilibrium 
conflguration of 'double-peaked' horns, where flux lost from 
the annulus is replaced by inflowing material. If the over- 
flow is sudden, the outermost of the double-peaked horns 
will decay rapidly as flux gets redistributed blue-ward and 
red-ward of the horns. Thus, if we observe a broad Fe K 
component with double-peaked horns that do not change 
over the viscous timescale, there is a dynamic equilibrium 
between pile-up behind and leakage from the dam. The en- 
ergy resolution of Astro-H is necessary to begin to detect the 
'double-peaked' horns corresponding to 'dammed' accretion 
disks. However, to systematically study such effects (rather 
than rely on serendipity), future proposed missions will re- 
quire much larger effective area than even that planned for 
IXO/Athena. 

4.3 Detecting a disk around the secondary black 
hole 

In our discussion of gaps and cavities above, we have ignored 
accretion onto the secondary. However (see i) I2.2|> . it is pos- 
sible that the secondary's accretion disk persists long after 
the inner disk has drained. An accretion disk within the Hill 
radius Ru = [q/iY^^R oi a, q — 0.03 secondary located at 
50rg in a cavity in the innermost AGN disk, has a radius 
of llrg in units of the primary black hole mass or 365rg,s, 
where rg^s corresponds to gravitational radii in units of the 
secondary black hole mass. The viscous timescale in the disk 
around the secondary is Ta,s = {l/as){Hs/rs)~'^l/ijJs where 
the parameters are as in eqn.([9]), but the s-subscript refers 
to the accretion disk around the secondary (which may have 
different properties to the main AGN disk around the pri- 
mary). Assuming ~ 0.01, //s/rs ~ 0.1, Ta,s ~ 10''rd,s, 
where rd,s is the dynamical timescale in the secondary disk, 
a disk extending to 365r-g,s around a g = 0.03 secondary 
orbiting a lO'^^^^M© supermassive primary black hole, will 
drain in ~ 2(200)yrs. Periodic dam overflows or continu- 
ous dam leakage could keep the disk persisting (see ^2.2|l 
around the secondary long after a cavity has formed. There 
could also be a small remaining inner disk around the pri- 
mary in this case, which will add flux to the red and blue 
wings of the Fe complex. Since motion will be about the cen- 
ter of mass of the merging binary, a small disk component 
around the primary will wobble, which may be detectable 
for large mass secondaries {q ^ 0.01) with the energy reso- 
lution of future missions such as Astro-H and IXO/Athena, 
however, we shall leave discussion of the effect of oscillations 
of the primary component to future work. In the following 
discussion, we consider only flux from simple accretion disks 
around the secondary and ignore the flux contribution from 
the gas streamers within the cavity. 
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Figure 11. The change in the observed broad FeK« Hne profile 
due to an accretion disk around the secondary black hole located 
at bOrg . The black curve corresponds to the Fe Ko emission from 
the main AGN disk (55 — lOOrg) plus a weak secondary broad 
component (10% of the intensity of the full disk profile) due to an 
accretion disk around the secondary black hole centered on the 
line centroid energy (6.40keV). The red curve shows the effect 
of shifting the centroid of the weak secondary component red- 
ward to 5.60keV. The blue curve shows the effect of shifting the 
centroid of the weak secondary component blue-ward to 7.05keV. 
The progression from red curve to blue curve occurs over half the 
orbital time of the secondary. Observationally, we expect a 'see- 
saw' oscillation between the blue and red wings of the line as the 
secondary accretion disk orbits the primary black hole inside the 
circumbinary disk cavity. All curves are binned at approximately 
the energy resolution {~ 7eV) expected for Astro-H. 



The secondary disk will add an additional component 
to the observed broad Fe Ka line as it orbits the primary. 
The observational effect will be most obvious at energies 
red- ward and blue- ward of the line component produced by 
the main AGN disk outside the cavity (e.g. the purple curve 
in Fig. [7]). A key difference between the ripple effect due 
to a disk gap, and the effect due to a mini-disk around the 
secondary, is that the former occurs on the AGN disk vis- 
cous timescale and the latter occurs on the much shorter 
orbital timescale. Thus, the effect consists of adding a sec- 
ondary, low intensity broad Fe Ka component that oscillates 
across the primary broad Fe Ka component on the orbital 
timescale. The secondary line component oscillates approx- 
imately between the red and blue horns of the primary line 
component over the orbital timescale. For a secondary lo- 
cated at 50rg from a IO^'^'Mq primary, one complete cycle 
of oscillation of the line centroid of the secondary compo- 
nent between the red and blue horns of the primary com- 
ponent will take ~ 2nRa/vs ~ lO'^'^-'s where Rs = 50rg is 
the secondary orbital radius and Vs is the secondary orbital 
velocity. 

In Fig. [11] we show the effect of adding an accretion 
disk around a secondary black hole located at 50rg in a 



cavity which truncates the AGN disk at 55rg. We assume 
that the circumsecondary disk properties are identical to the 
circumbinary disk properties (emissivity, disk inclination) 
except for the cavity. We also assume that the secondary 
broad line flux is 10% of the primary flux without a cavity 
(the black curve in Fig. [S]). In Fig. Illl the black solid curve 
corresponds to the superposition of the broad line proflle 
from a AGN disk with cavity (purple curve in Fig. [7| , with 
a secondary broad line component centered on the line cen- 
troid energy (6.40keV). The red-curve in Fig. Illl shows the 
AGN disk/cavity profile with the secondary disk line profile 
superimposed at the most red-shifted part of its orbit (line 
centroid is ~ 5.6keV). The blue curve in Fig. Illl shows the 
AGN disk/cavity line profile, with the secondary disk line 
profile superimposed at the most blue-shifted part of its or- 
bit (line centroid is ~ 7.05keV). Thus, Fig. [11] shows the 
effect of one half of the orbit, as the secondary component 
centroid moves between the red horn to the blue horn of 
the primary component. The key observable for detecting 
disks around secondary, satellite black holes is therefore a 
periodic see-saw oscillation between the red and blue curves 
as in Fig. Illl If there were a weak circumprimary disk (say 
^ lOrg), it would add a very broad component to the red 
and blue wings of the Fe K line complex. Observationally 
this would be equivalent to adding a component to the red 
and blue wings that will vary on the timescale of the pri- 
mary wobble around the binary center-of-mass. The wobble 
of the primary component depends on the mass ratio (q) of 
primary to secondary and may be detectable with the en- 
ergy resolution of IXO/Athena in very long exposures. Here 
we shall focus on oscillations due to the secondary disk and 
we shall leave discussion of the effect of oscillations of the 
primary component to future work. 

In Fig. [12] we show the same effect, this time for a sec- 
ondary located at 30rg inside a cavity that truncates the 
disk at 55rg. In this case, we see a much more pronounced 
oscillation between the red and blue horns of the total ob- 
served broad line component. Here, the line centroid energy 
of the secondary oscillates between 5.2keV and 7.3keV. This 
oscillation can occur on the very short secondary orbital 
timescales (few-lOs of ks in AGN with a primary 10^~^ Mq) 
and could easily be detected during an extended observation 
with Astro-H. Thus, the deeper the secondary accretion disk 
in the cavity, the easier it is to observe the osci llation effect. 
Fig- dUlooks quite different from the results of lSesana et al.l 
(i2012l' ) (their Fig. 10), where there are two LAOR model disk 
components, corresponding to a W^Mq total-mass binary, 
with the secondary blue-shifted by lO^kms"^ (correspond- 
ing to a separation of ~ lO^rg), much larger than the sep- 
arations discussed here. The apparent differe nce is mostly 
driven by the lack of horns in the profiles in ISesana et al.l 
(|2012| ) which is due to thei r choice of steeper emissiv ity pro- 
file r"^ (see e.g. Fig. 10 in ([Reynolds fc Nowakllioosh for the 
impact of this choice on broad line profile) . 

In Fig. [TH] we simulat«lf] a 350ks observation with 
XMM-Newton EPIC-pn of an AGN at 2 = 0.01 with 
a 2-lOkeV photon fiux of 3.5 x 10""photcm~^s"\ We 
ignore the narrow line component in the model in or- 
der to illustrate the difference in best-fits to observa- 



http:/ /heasarc. gsfc.nasa.gov/cgi-bin/webspec 
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Figure 12. The change in the observed broad FeK« Hne profile 
due to an accretion disk around the secondary black hole as in 
Fig. 111! except the secondary is located closer in, at 30rg. The 
black curve corresponds to the Fe Ka emission from the main 
AGN disk (55 — lOOrg) plus a weak secondary broad component 
(10% of the intensity of the full disk profile) due to an accretion 
disk around the secondary black hole centered on the line centroid 
energy (6.40keV). The red curve shows the effect of shifting the 
centroid of the weak secondary component red-ward to 5.2keV. 
The blue curve shows the effect of shifting the centroid of the weak 
secondary component blue-ward to 7.3keV. The progression from 
red curve to blue curve occurs over half the orbital time of the sec- 
ondary. Observationally, we expect a 'see-saw' oscillation between 
the blue and red wings of the line as the secondary accretion disk 
orbits the primary. All curves are binned at approximately the 
energy resolution (~ 7eV) expected for Astro-H. 



tions; in practice this model component would need to 
be fit to the data too. The XSPEC toy model was sim- 
ply ZPHABS(zPOWERLW+DlSKLlNE-f diskline) where the ab- 
sorbing Hydrogen column was lO^^cm"^ and the power- 
law index was —1.8. The primary diskline was centered at 
6.40keV and due to a truncated inner disk from 50-lOOrg, 
with 9 = 60° and emissivity k = —2.5 as above, and a line 
normalization of 10~^ with respect to the continuum. The 
secondary diskline had similar parameters to the primary 
diskline except it ran from 6-100rg,a, where rg,s represents 
gravitational radii of the secondary. We set the normaliza- 
tion of the secondary diskline to ~ 0.2 of the primary to 
illustrate the difference between the best fits to two simu- 
lated observations. The simulated data in in Fig. [13] corre- 
sponds to line centroid of the secondary at 7.3keV and the 
blue solid line is the best model fit to the data. The red 
solid line in Fig. [13] corresponds to the best model fit when 
the line centroid of the secondary lies at 5.2keV, although 
we do not show the simulated data in the latter case for 
ease of presentation. Flux differences in the red wing be- 
tween the two models in Fig. [TJ] differ by a few percent and 
the key observational test in this case is that the best-fit 
model must oscillate between the blue and red curves on a 
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Figure 13. A simulated 350ks observation with XMM-Newton 
EPIC PN of an AGN at ^ = 0.01, with a 2-lOkeV count- 
rate of 3.5 X 10~^^photcm~^s~^, which is best fit with a 
two component diskline model (assuming the narrow FeKa 
component has been subtracted). We used a simple XSPEC 
model zphabs(zpowerlw+disklineH-diskline) where the ab- 
sorbing Hydrogen column was lO^^cm"^ and the powerlaw in- 
dex was —1.8. The primary diskline has line centroid at 6.40keV, 
emissivity k = —2.5, d = 60° and originates from 50 — lOOrg in a 
truncated inner disk, with line normalization set to 10~^ of the 
continuum. We ignore the narrow Fe Ka line component for clar- 
ity.The secondary diskline has similar parameters to the primary, 
but relative normalization of 0.2 in order to highlight the sense 
in which the oscillation can change the best-fit. Black data points 
correspond to simulated data where the secondary line centroid 
lies at 7.3keV and the blue solid line corresponds to the best-fit 
model to the data. The red solid line corresponds to the best 
model fit when the line centroid of the secondary lies instead at 
5.2koV (simulated data not shown for clarity). 



regular orbital timescale. As we can see, even in two long 
simulated exposures with XMM-Newton and a relatively in- 
tense secondary disk component, it is difficult to disentangle 
the secondary disk component, even if we bin up the data. 
However, repeated observations of bright AGN can build up 
statistical significance in the red and blue wings, allowing 
us to track the oscillation of the secondary over time. Thus, 
observations of AGN in the Fe K band with XMM-Newton 
EPIC-PN should be searched for serendipitous oscillations. 
A systematic study of models of secondary disks in AGN 
could be carried out with LOFT since the proposed effec- 
tive area (lOm'^) is an order of magnitude larger than XMM- 
Newton and the spectral resolution (~ 0.25keV) is sufficient 
for coarse binning of the red and blue wings of the broad 
Fe Ka line. For disentangling competing effects, ideally a 
Large Area High Resolution (LAHR) telescope should be 
used, with effective area greater than that of LOFT together 
with the proposed spectral resolution of 1X0/ Athena. 
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5 MULTIMESSENGER ASTRONOMY WITH 
GRAVITATIONAL WAVES AND THE 
RIPPLE EFFECT 

Observations of a ripple effect or periodic oscillations in the 
broad FeKa line will be an excellent leading indicator of 
gravitational radiation which can be detected directly in- 
dependently. Depending on the orbital period of the bi- 
nary, the GWs emitted by a supermassive black hole bi- 
nary are detectable with pulsar timin£| or a space-based 
laser-interferometric instrument like LISA[f| or NGC|f|. The 
GW measurement can constrain the source parameters com- 
pletely independently of the electromagnetic signal with dif- 
ferent systematics. Furthermore, the GW measurement may 
independently indicate the presence of ambient gas around 
the binary as discussed below. A serendipitous binary detec- 
tion in one channel can be used to do a tar geted search in the 
other for verification (|Kocsis et al.|[200^ ). The comparison 
of these measurements can provide a consistency check for 
the central mass, inclination, and orbital radius. The combi- 
nation of the electromagnetic and GW measurements may 
allow us to constrain the geometry of the accretion disk and 
the binary. In this section we discuss the gravitational wave 
signatures that may be coincident with the electromagnetic 
signatures. 

The GW strain amplitude generated by a circulaiQ bi- 
nary of separation r, distance D, inclination t relative to the 
line of sight is 



h+{t) 



= 2 



mi 7712 

rD 



(1 + cos^ l) sin20(t) 
2 cos L COS 2(j){t) 



(10) 



where (j){t) is the orbital phase along the orbit. Ifere, the 
separation r changes in a characteristic way due to GW 
losses and gas effects, which causes changes in the orbital fre- 
quency, d<j)/dt. The corresponding phase shift can be used to 
measure the binary paramete rs including the separation , ec- 
centricity, masses, and spins (|Cutler fc FlanaganllT993 ) and 
the torques exerted on the binary by the accretion dis k 
(|Kocsis et al.ll201ll : lYunes et al.ll201ll : iHavasaki et"aDl2012h . 
The distance to the source and the binary inclination then 
follows from the amplitudes of the two GW polarizations 
h+ and hx- The sky position of the source can be inferred 
by measuring the GW strain coincidentally over multiple 
baselines. 



5.1 PTA detections 

Pulsar timing arrays (PTAs) can measure GWs of periods 
weeks to years (or frequencies between nHz to mHz) by de- 
tecting correlated variations in the time-of-arrivals (TOAs) 
of pulses in millisecond pulsars in the Galaxy. At these pe- 
riods the binaries are relatively far from merger, so only 
four quantities are measurable for most of the sources: the 
GW amplitudes of the two polarizations (see Eq. I10[) . the 
orbital period, and the eccentricity, but not the inspiral rate 



^ 'http : //www . ipta4gw . org/ ' 
^ http://lisa.nasa.gov/ 
^ |http: / /elisa-ngo . org/ 

''' The GW waveform for general eccentric orbits is sim- 
ilar, but additional ter ms involving (p appear in that case 
jPeters fc Mathews|[l963l ). 



or <f>. Ho wever, the inspi ral rate may be measurable for some 
sources (|Lee et al.ll2"oil] ). In the case of PTA detections, we 
expect to dete ct the binary when it is close to decoupling 
from the disk jSochora et al.llioill 'l and therefore the signa- 
ture of a gap or cavity in the broad Fe Ka line caused by 
disk drainage (outlined above) should be clearest. We expect 
to detect a few such binaries with PTAs out to z < 1. 

These observations are mostly sensitive to the total 
stochastic background generated by many supermassive 
black hole binaries acr oss the Uni verse at low frequencies 
(.Phinnev .2001: Sesan a et al.ll2008l l. However the GWs of 
massive nearby binaries (e.g. mi ~ m2 > 10^ Mq within 
z < 1-2) are expected to rise well ab ove the background 
level and may be indiv idually resolved (jSesana et al.l l 20091 : 
iKocsis fc SesanallioTll ). The individually detectable bina- 
ries are expected to be orbiting at separations 100 — 200rg, 
where gas-driven mi gration is subdominant relative to gravi- 
tational wave losses ( Kocsis fc Sesanallionl ). These GW ob- 
servations would be able to localize such sources on the 
sky by comparing the TOAs of several pulsars in differ- 
ent directions. The sky localization and inclinatio n measure- 
ment accuracy are of order degrees llSesana fc Ve cchio 20 13; 
iLee et al.|[2oTll : iBabak fc Sesanall2012l : |Petiteau et al.ii201 . 

The individually resolvable massive binaries have a typ- 
ically comparable mass-ratio {q ~ 1). These binaries are 
expected to open a large cavity in the disk. The FeKa pro- 
files may show a double structure if the non-axisymmetric 
streams can penetrate the gap a nd create an accretion disk 
around the binary components (jSesana et al.li20l2 '). These 
bin aries may also exhibit periodic electromagnetic variabi l- 
ity (|MacFadven fc Milosavlievidl2008l : iHaiman et al.ll2009l '). 
With large cavities expected to be excavated in the inner 
disk, we expect to simultaneously detect broad Fe Ka pro- 
files in these AGN similar to the individual curves shown 
in § m If we are lucky enough to detect a massive binary 
via PTA as the cavity is forming, we will simultaneously ob- 
serve the variable broad Fe Ka profiles, together with ripple 
effects and oscillations outlined in § |4] 

Since the inspiral rate is typically not available for these 
sources, one cannot directly measure the effects of gas for 
individual binaries from PTA observations. Gas driven mi- 
gration may leave an imprint on the stochastic GW back- 
ground by reducing the number of binaries in the Universe 
at large separations l|Kocsis fc Sesanall201ll '). For individual 
binaries, the effects of an accretion disk may be inferred 
from the GW signal indirectly by measuring a large orbital 
eccentricity and the absence of spin-orbit precession. The 
angular momentum exchange with an accretion disk leads 
to the excitation of orbital eccentricity for widely separated 
binaries if a cavity forms in the disk, while the eccentricity 
is predicted to quickly decrease in the absence of an accre- 
tion disk or during the later stages of the inspiral due to 
GW emission. Ho wever, the eccentricity is also increased b y 
scattering of stars l|Amaro-Seoane et al.ll2009l : ISesanair2010l ) , 
albeit to a smaller level. The eccentr icity may be detected 
with PTAs (|Enoki fc Nagashimall2007l 'l. SMBH binaries with 
a large orbital eccentricity may be indicative of an accre- 
tion disk. Secondly, spin-orbit precession may also be de- 
tect able if the spins ar e not aligned with the orbital plane 
( Mi ngareUi et al.1 120121 ') . Such a configuration would argue 
against the presenc e of a gas disk, which tends to align the 
binary spins (Bogdanovic. Revnolds fc Milled I2OO7I ). 
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5.2 LISA/NGO detections 

LISA/NGO will be sensitive to much smaller periods be- 
tween 10 hours to seconds, which corresponds to frequencies 
from 3 X 10"^ to IHz and SMBH masses of M ~ 10^- 
10^ Mq at separations from ~ 100 Vg down to merger. In 
this frequency range, the inspiral of SMBH binaries is much 
more rapid and the astrophysical stochastic GW background 
is muc h more quiet, dominate d by Galactic white dwarf bi- 
naries (jNelemans et al.ll200lh . In this case, for a massive 
binary very close to merger, if an annulus is maintained in 
the disk, we expect to simultaneously observe a broad Fe Ka 
component in the late stage of the ripple effect (e.g. the red 
curve in Fig. [l]). If instead, the binary merger is occurring 
within a cavity, we expect a strongly dammed cavity inner 
wall, and a simultaneous Fe Ka broad line profile similar to 
the green curve in Fig. [TO] 

The GW signal of an IMBH or an SMBH spiraling 
into a SMBH is much stronger than the background level, 
and may be detected to high significance to high redshifts 
(z ^ 10 for comparable m ass ratio binaries, Hughes 2002; 
lAmaro-Seoane et al.ll201^ ). A sensitive measurement of the 
inspiral waveform for the two GW polarizations can be 
used to infer the binary masses and orbital parameters to 
a relative accuracy better than 10~^ for a signal to noise 
of 30 (jfiarack & Cutler 20o3). The source localization ac- 
curacy with LISA/NGO is typically a few d egrees sev- 
eral months before merger (|Kocsis et al.l [ioO^ l. This im- 
proves to a few tens of arcminutes for comparable-mass 
spinning binaries approaching the I SCO at z ^ 1 if the 
spins are not aligned w ith the disk (|Lang fc Hughe jl2008l : 
iMcWilliams et ani201lt). If the disk plane and t he binary 
spin align (|Bogdanovic. Revnolds fc Milled [2OO7I ) . there is 
no spin-orbit precession feature in the GW signal, and the 
binary localizatio n becomes worse by a factor of a few 
jLang et al.ll201ll ). This error volume however is still suf- 
ficiently sma ll to allow a uniqu e identification of an AGN 
counterpart (|Kocsis et akllioOB ). 

Due to the high signal to noise ratio and the large num- 
ber of orbital cycles observable with LISA/NGO (A^cyc ^ 
10^), the imprints of gas effects may be detected directly 
by me asuring the corresp o nding perturbatio ns of the GW 
phase l|Kocsis et al.ll201ll : lYunes et all 120111 ). The largest 
perturbation is due to the spiral density waves excited by 
the binary which removes angular momentum from the bi- 
nary analogous to planetary migration. Further perturba- 
tions arise due to the hydrodynamic drag on the secondary 
and the mass increase of the secondary due to accretion. 
These effects are more significant for smaller secondary 
masses. For an extreme mass ratio binary, mi/m2 < 10~*, 
the GW phase shift due to the spiral waves are expected 
to be measurable w i th LISA during th e last year of inspiral 
(jKocsis et al.ll201ll : FVunes et aLlbOllT ). The GW measure- 
ments may be directly sensitive to the presence of a gap 
since a gap affects the gaseous torques acting on the binary, 
and hence, the GW phase. 

For comparable mass-ratio binaries, gas effects are much 
less significant, only detectable by second ge neration space- 
based instruments like Decigo or BBO (jHavasaki et al.l 
I2OI2I ). In this case, the merger may take place in a cen- 
tral empty cavity in the disk where the gaseous torques 
are negligible. However, eccentricity of order e = 10~3 



to 10^"^ in LISA/NGO detections may be an indirect sig- 
nature of an earlier coupling with a massive circumbi- 
nary disk tArmitage fc Na taraian 2005; Roedig ct al. 2011|). 
LISA/NGO will allow to measure the eccentr i city to a much 
higher accuracy 10-* ijBarack fc Cutleill2004l : [Mikoczi et al.l 



6 CONCLUSIONS AND FUTURE WORK 

If a gap-opening, migrating black hole ends up in the inner- 
most regions of an AGN disk, analogous to a gap-opening 
'hot Jupiter' in a protoplanetary disk, a unique and pre- 
dictable pattern of variability appears in the broad com- 
ponent of the Fe Ka line. The ripple effect and oscilla- 
tions in the Fe Ka line outlined in this paper are poten- 
tially detectable in long exposures with future missions such 
as Astro-H, IXO/Athena and LOFT. Oscillations may be 
detectable in long, repeated exposures with XMM-Newton 
EPIC-PN, such as may be found for popular sources in the 
archive. Here we have shown the range of ripple effects ob- 
servable as a function of disk and merger properties, includ- 
ing inclination to observers' sightline, gap-width (or cavity 
size), disk emissivity profile, damming of the accretion flow 
and a mini-disk around the secondary black hole. 

Detection of a ripple effect or periodic oscillations in 
the broad component of Fe Ka from an AGN will provide 
advance warning of gravitational waves due to an impending 
merger in this AGN. For example, an observation of a broad 
FcKq line ripple profile given by the dashed red-curve in 
Fig.[l]from a supermassive black hole of mass ~ IO^Mq pre- 
dicts a merger event in this source within a year. Once grav- 
itational waves consistent with a binary black hole merger 
are detected, an archival search for a ripple effect or oscilla- 
tions in a broad FeKa line will help localize the gravitational 
wave detection. Departures from the predicted ripple effect 
in the final stages of merger will allow us to test the pre- 
dictions of strong gravity independent of the detection of 
gravitational radiation. 
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